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AN  EXTENSION  THEOREM  FOR  FINITE  ELEMENT  SPACES 
WITH  THREE  APPLICATIONS 

Olof  B.  Widlund 
Courant  Institute  of  Mathematical  Sciences 
251  Mercer  Street 
New  York,  NY  10012   USA 

SUmiARY 

In  this  paper,  an  extension  theorem,  siiriilar   to  well-known 
results  for  Sobolev  spaces,  is  established  for  general  conform- 
ing finite  element  spaces.  Results  of  this  type  provide  central 
tools  for  the  theory  of  iterative  substructuring  (domain 
decom.position)  and  capacitance  m.atrix  m;ethods.   The  main  theorem, 
has  previously  been  established,   using  other  techniques,  for 
certain  Lagrangian  finite  elem.ents  in  the  plane.   Several  appli- 
cations are  discussed  in  which  previous  results  on  iterative 
substructuring  and  capacitance  m.atrices  are  extended  to  three 
dimensions,   non-Lagrangian  finite  elements   and  higher  order    , 
elliptic  problems. 


1.   INTRODUCTION 

Recently  there  has  been  a  considerable  interest  in  the 
developm.ent  and  study  of  iterative  substructuring  methods  for 
elliptic  finite  elem.ent  problem.s;  see  Bjizirstad  and  Widlund 
[6,7],  Bram.ble,  Pasciak  and  Schatz  [8,9],  Chan  [11]  ,  Chan 
and  Reasco  [12],  Dihn,  Glowinski   and  Periaux  [16]  ,    Dryj a 
[17,18],  Dryja   and  Proskurov/ski  [19,20],  Golub  and  Mayers  [21], 
Keyes  and  Gropp  [22]   and  l^idlund  [29].   Early  work  in  this 
area  is  described  in  Concus,  Golub  and  O'Leary  [15],  while  the 
paper  by  Dryja   [17]   appears  to  be  the  first  in  v/hich  an 
optimal   iterative  method  is  described  and  analyzed.   Much  of 
this  work  has  been  for  finite  difference  approximations  of 
second  order  elliptic  problems,  often  on  relatively  special 
plane  regions.   It  is  the  purpose  of  this  paper  to  go  beyond 
these  cases  to  more  difficult  and  important  finite  element 
problems . 

When  using  these  algorithms,  also  known  as  domain  decomp- 
osition methods,   the  discretized  elliptic  problem  is  parti- 
tioned into  subproblems  which  correspond  to  non-overlapping 
subsets  of  the  region.   The  subproblems  are  then  solved 
separately,  and  repeatedly,  while  the  interaction  between  the 
subregions  is  handled  by  a  conjugate  gradient  or  other  suitable 
iterative  method.   These  iterative  methods  provide  interesting 
alternatives   to  the  standard  industrial  finite  element  practice 
in  which  not  only  the  stiffness  m.atrices  corresponding  to  the 
finite  element  m.odels  for  the  substructures  but  also  the 


matrices  which  represent  the  interaction  betv/een  the  different 
parts  are  fully  assembled  and  factored  into  their  Choleski 
triangular  factors;  see  Bell,  Hatlestad,   Hausteen  and  Ar  Idsen 
[2  ]   and  Przemieniecki  [25].   For  a  discussion  of  different 
algoritnms  etc.  and  a  general  finite  element  framework,  see 
Bj0rstad  and  Widlund  [  7 ] . 

The  fastest  among  these  methods  offers  considerable 
advantages  even  on  sequential  computers.   They  also  shov; 
particular  promise  for  parallel  computing  since  the  last  stage 
of  a  block  Gaussian  factorization  is  likely  to  lend  itself  less 
well  to  parallel  architectures  than  the  earlier  stages  v/here 
work  can  be  carried  out  for  the  different  substructres  without 
any  need  of  communication  between  them.   Work  on  the  actual 
parallel  implementation  of  these  methods  is  beginning;  see 
Keyes  and  Gropp  [22]  for  results  on  an  Intel  Hypercube. 
Similarly,  variants  of  methods  for  general  sparse  linear  systems 
of  algebraic  equations,  such  as  nested  dissection  algorithms, 
could  be  developed  in  which  the  factorization  would  be  stopped 
at  a  suitable  stage.   The  remiaining  reduced  system,  correspond- 
ing to  the  variables  not  eliminated,  would  then  be  solved  by 
an  iterative  method. 

The  rate  of  convergence  of  iterative  methods  of  this  kind 
can  of  course  be  studied  by  conducting  systematic  numerical 
experiments.      For  general  sparse  or  band   matrices   it   is 
unlikely  that  strong  general  convergence  results  for  powerful 
preconditioners  Kil_  be  developed  for  general  system.s  of  equa- 
tions.  When  a  problem  arises  from  the  discretization  of  an 
elliptic  problem,   additional  tools  from,  mathematical  analysis 
are  available  and  the  possibility  of  a  systematic  developm^ent 
of  good  preconditioners  can  become  a  reality.   Thus  for  an 
elliptic  problem   with  constant  coefficients  on  a  uniform  mesh 
and  a  region  which  is  a  union  of  a  few  rectangles,  Fourier 
series  can  be  used;  see  e.g.  section  4  of  Bj^rstad  and  Widlund 
[  7  ]  .   Similar  techniques  are  used  in  Bj0rstad  [4,5]  to  study 
a  biharmonic  problem  and  in  Chan  [10]   to  develop  and  analyze 
new  fast  methods  for  queuing  networks.   However,  a  system.atic 
theory  for  elliptic  problems  on  general  regions  requires  the 
development  of  finite  element  analogues  of  certain  regularity 
results  for  inhomogeneous  elliptic  boundary  value  problems, 
see  Bj(Z$rstad  and  Widlund  [  6  /  7  ]  and  Bramble,  Pasciak  and 
Schatz  [8,9].   The  technical  aspects  of  this  work  is  compli- 
cated by  the  fact  that  the  boundaries  of  the  subregions 
necessarily  have  corners. 

As  in  the  continuous  case,  we  can  reduce  an  elliptic 
problemi  with  nonhomogeneous  boundary  data  to  one  which  is 
homogeneous  by  using  an  extension  theorem.;  see  e.g.  Stein  [21]. 
In  recent  papers,  see  Bj^rstad  and  Widlund  [  7  ]   and  Bram.ble, 
Pasciak  and  Schatz  [  8  ] ,   an  extension  theorem  has  been 
developed  for  conforming  Lagrangian  finite  elements  in  the 
plane.   This  proof,  which  essentially  has  been  knovm  at  least 
since  1980,  uses  a  variant  of  the  Bramble-Hilbert  lemma;  see 
Ciarlet  [13],   and  a  regularity  theorem  for  elliptic  differen- 
tial equations.   It  appears  difficult  to  extend  this  proof  to 
non-Lagrangian  finite  elements,  three  dimensions  and  higher 
order  equations  since  it  relies  on  the  boundedness,  in  a 


Sobolev  space,  of  the  finite  element  interpolation  operator. 
The  choice  of  this  Sobolev  space  is  limited  by  the  lack  of 
regularity  of  the  solutions  of  elliptic  problems  on  regions 
with  corners.   In  this  paper,  we  therefore  proceed  differently, 
inspired  by  an  idea  used  by  Astrakhantsev   [  1  ]   to  establish 
that  a  capacitance  matrix  method  for  a  rather  special  varia- 
tional difference  scheme  is  optim.al  for  the  Neumann  boundary 
condition. 

We  note  that  the  extension  theorems  are  related  to  questions 
concerning  the  approximation  order  of  finite  element  spaces  for 
functions  which  are  not  sufficiently  smooth.   These  issues  were 
discussed  early  by  Strang  [28]  for  cases  without    boundaries. 
This  work  was  extended  to  some  extent  by  Clement  [14]   and  later 
quite  systematically  by  Bernardi  [  3  ]  .   VJhile  Bernardi  carefully 
includes  boundary  conditions  in  her  work,  her  work  does  not 
include  Hermitian  finite  elements   and  it  is  therefore  not 
directly  applicable  to  all  the  problems  at  hand. 

In  section  3,   v/e  discuss  three  applications.   Ke  thus 
extend  Astrakhantsev' s   result   on  capacitance  matrix  methods 
for  Neum.ann  problem.s  to  a  much  more  general  family  of  finite 
elements.   While,  at  least  in  the  West,  iterative  substructuring 
methods  are  attracting  much  more  attention  than  capacitance 
m.atrix  algorithms   the  tv^o  famalies  of  methods  are  nevertheless  - 
closely  related.   We  also  give  two  new  results  on  domain 
decomposition  in  this  last  section. 


2.   THE  EXTENSION  THEOREM 

Let  Q.   be  an  open,  bounded  region  in  R  ,  d  >_  2 ,  with  a 
piecewise  smooth,  uniformly  Lipschitz  continuous  boundary  T. 
It  is  well  known,  cf.  e.g.  Stein  [27],  Chapter  VI,  that  there 
exists  a  linear  operator  ft  extendina  functions  on  2  to  functions 
on  R^  with  the  properties 


i) 


ftul^  =  "' 


i.e.  £  is  an  extension  operator, 

i.e.  £  maps  H^  {9.)    continuously  into  h'^CR  ) 
integer. 


(2.1) 


m  is  a  nonneaative 


The  constant  C(fi)  depends  on  d,  m  and  the  Lipschitz 


constant  of  the  region  only.   The  Sobolev  space    H  (o)  is  the 
subspace  of  L_ (o)  such  that 
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and  I  a  I  =  ^  a . . 


We  will  develop  an  analogue  of  this  extension  theorem  for 
finite    element  spaces.   Let   t^  be  a  triangulation  of  Q  into 
elem.ents.   For  convenience  we  assume  that  each  element  K  is  a 
simplex  with  sufficiently  smooth  faces;  our  theory  could  equally 
well  be  developed  for  quadrilaterals.   The  simplices  can  be 
curved.   As  usual  we  assume  that  the  simplices  are  properly 
joined  which  means  that  the  intersection  of  any  two  of  them, 
is  empty  or  an  entire  face  of  dimension  d-1  or  less.   By  hj^  v/e 
denote  the  diameter  of  K  and  by   pj^  the  diam.eter  of  the  largest 
ball  contained  in  the  element.    Following  Bernardi  [ 3  ] ,  we 
assume  that  each  element  is  the  image  of  a  straight  d-simplex  K 
under  a  C^  map  F-^      of  the  formi 

P  -   -p      +  $ 
K     K     K 

where   F„  is  an  invertible  affine  mapping 

K 
F^:    X  -  ij^x  +  h^ 

and  $„  satisfies 

K 

c^   =   sup  I  \D<i^{x)B^    II  <  c  <  1 

with  c  uniformly  bounded  away  from  1.   VJe  note  that  K  =  F„(K) 
is  a  straight  d-sim.plex  and  that  we  can^arrange  it  so  that 
K  and  K   have  the  same  vertices,  i.e.   K  is  a  straight  approxi- 
mation of  K.   Isoparametric  elements  as  well  as  those  based  on 
an  exact  triangulation,  cf.  Bernardi  [  3  ],  Lenoir  [23]   and 
Scott  [26],  can  easily  be  accommodated  in  this^f ramework.   The 
triangulation  is  assumed  to  be  regular,  i.e.   ^k^'^K  —  °  '  ^^^^s 
o  is  uniformly  bounded  for  all  h.   These  assumptions  exclude 
degenerate,  very  flat  elements  and  also  lead   to  a  bound  on  the 
number  of  simplices  that  have  a  particular  vertex  in  common, 
but  still  allow  us  selectively  to  refine  the  triangulation  to 
improve  the  accuracy  where  the  solution  is  less  smooth.   In 
practice  a  simtplex  will  be  straight,  unless  at  least  tv.'o  of 
its  vertices  fall   on  the  boundary  or  an  interface  between 
subregions  modeled   separately  e.g.   because  of  different 
material  properties.   We  will  not  discuss  the  effects  of  approxi- 
mating Q  by  a  union  of  straight  of  isoparametric  elements  but 
always  assume  that  the  triangularions  of   fi  are  exact. 

We  will  consider  general,  conforming  finite  elements,  i.e. 
the  approximating  space  V^  c  v  where  V  is  the  linear    space 
appropriate  for  the  elliptic  variational  problem  to  be  consid- 
ered; see  Ciarlet  [13].   A  finite  element  is  defined  on  the 
element  level  by  a  triple  (K,Pj^,Zj^)  where  K  is  the  simplex  and 
Pj.  a  space  of  functions.  The  dimension  of  Pj^  is  bounded  uni- 
formly.  Normally  Pj.  is  the  image  under  the  map  Fj,  of  a  space 
of  polynomials.   We  note  that  these  functions  can  be  vectored 
values.   No  essential  difficulties  are  introduced  by  such  an 
assumption  nor  by  considering  curved  rather  than    straight 
simplices.  Z-^   is  a  set  of  linear  functionals  defined  on  smooth 
functions  and  on  Pj^  such  that  the  related  interpolation  problem 
on  K  is  uniquely  solvable  in  P-^. 

From  the  set   ^    Z^,  a  m.aximal  system  of  linearly  inde- 


pendent   linear   functionals  on   C  (il)    is  extracted.   Basis 
functions  for  the  finite  element  space   v"  can  then  be  con-  . 

structed  as  e.g.  in  Ciarlet  [13],  creating  {(t-}-_-,   and  {li'}-_-| 

which  are  a  dual  pair  of  bases  for  the  space  V"  and  the  space 
of  linear  functionals,  respectively   i.e. 

^i^*j^  "  *ij  '    ^   1   ^   1\    '         ^  1   ^    1\    ' 
For  1  <_  i  <_  Nj^  ,  we  define 

A.  =   u    {K;  supp  y .  overlaps  K}  . 
^    KGxh  ^ 

Under  the  assumptions  given  above  one  can  show  straightforwardly 
that  the  number  of  elements  that  form  A_t  is  uniformly  bounded 
and  that  the  diameters  of  any  elements  K  and  K' ,  contained  in 
the  same  Aj_  satisfy 

where  C   depends  on  c  =  max  Cj^  and  o  -   max  hj^/pj^  only.   The  sup- 
port of  (|:j_  is  contained  in  Aj_  . 

In  our  proof  we  assume,  following  Strang  [28],  that  the 
basis  functions  (+)j_  (x)  are  uniform  of  order  m.,  i.e.  there  exists  - 
constants   Cg  ,  independent  of  h  and  i,  such  that 

max  |(|-)^d).(x)|  <  c   h  ^    ,   s  <  m  .  (2.2) 

xSA.   ^"^    ^      "   ^ 

1 

|a!=s 

We  note  that  dj^  often  is  the  degree  of  a  derivative  associ- 
ated with  the  degree  of  freedom  in  question. 

Assuming  that  CT. ,  the  complement  of  r2,has  been  triangulated 
in  an  equally  benign  way  and  that  v'^C:^)  has  been  extended  to 
V^(R^)  c  H^CR*^),  we  are  ready  to  describe  how  to  find  an  exten- 
sion S^uj^  G  vh(Rd)   of  a  given  elem.ent   u^  ^  v"(.T)   such  that 

M^\i!  m   d  1  Cl|u  II       .  (2.3) 

In  a  first  step,  we  use  the  original  extension  theorem  and 
the  fact  that  the  space  is  conforming,  i.e.  u,  ^  H^(n),  to  find 
&u^   e  h"''(R'^)   such  that 

I  |Su,  I  I      ,   <  C(r2)  I  |u,  I  I       .  (2.4) 

'    "^   h"^(R^)  -         ^   H"^(fi) 

However  &u^  ^   V  (R  ). 

In  a  second  step,  we  follow  Strang  125]  and  smooth  Su  and 
interpolate  to  find  an  element  w,  e  V  (R  )  such  that 

l|w,  I  1     .   <  C||£u,  I  I     .   <  C  (n)C  Iju,  I  I       .    (2.5) 

"  h  h^(r'^)  -     ^  H^(R^)  -       ^  H"(a) 

However  w,  |  ^  7^  u,  and  w,  is  therefore  not  an  extension  of  u,  . 


We 


note  that  Strang   [28]  shows  that 


h,   ^     ^,  ^  h 


l"h 


-    f^u,     T     <  Ch„       }  w,  -£u,  .      (2,6, 


where  the  suir,  is  over  a  fixed  number  of  neighboring  eleirients. 

In  a  third  step,  we  patch  w-^    and  u^    together  to  create  £  u,  . 
We  interpolate  us_ing  all  the  original_paraineters  obtained  from 
Uy^  in  Q,   and  in  CQ,  the  cpmplem.ent  of  Q,  those  associated  with  w.  . 
The  resulting  function  S  u,  e  v^  (R'^)  ,   For  any  element  K  c  si 
the  original  values  of  u^  are  recovered.   We  also  note  that 
similarly,  in  CQ ,    S  u^  differs  from  w^  only  on  elements  which 
have  at  least  one  vertex  on  F.   There  remains  to  establish  (2.3). 

By  the  triangle  inequality 

I  K-'^u,  I  I       <  ||8^u,-w,  II        +  l|w,  II 

By  (2.1)  and  (2.5),  |  |  w,  |  |  „,       can  be  estimated  bv  |  |  u,  |  | 

'  '  hi  Ih"^(CQ)  -  '  '  h'^m^^j 

Consider  an  element  K'  c  CQ,    with  at  least  one  vertex 
on   r  for  which  S'^uj^^-Wj^  is  not  identically  zero.   The  parameters 
associated  with  the  local  interpolation  problem  are  equal  to 
zero  except  for  those  shared  with  at  least  one  elemient  K  C  f^. 
By  using  (2,2)  one  can  establish  that  "7 

||£^u,-w,  ||2^      <  ch'^'^l  |£^u,  -w,  I  |^_      .  (2.7) 

h   hi  '^'^(j,.)  -    K   i  I    h   hi  Il2(k') 

Under  the  regularity  assumptions  introduced  above  and  (2.2), 
we  can  estimate  this  L2  norm  by  h„  '  xthe  I2    norm  of  the  vector 
of  parameters  h*-^iy  •  (£"u,  -wj^)  .  When  considering  the  correspond- 
ing norms  of  S^u,  -  w,   on  the  neighboring  element  K  C  fi ,   we 
see  that 

I  I  £  u,  -w,  I  I  T      <  C  I  I  £  u,  -w,  I  I  T     =  C  I  I  £  u,  -w,  |  |  ^ 

The  crucial  observation  is  that  this  inequality  holds  whatever 
the  values  of  the  additional  parameters  associated    with   K 
happen  to  be. 

The  proof  is  concluded  by  adding  over  all  the  relevant 
elements  of  CQ   and  using  inequalities  (2.6)  and  (2.7).   An 
exact  analogue  of  the  extension  theorem  given  in  the  beginning 
of  this  section  has  thus  been  established  for  conform.ing  finite 
element  spaces  defined  on  sufficiently  regular  triangulations . 
We  note  t^at  £   is  a  linear  operator  and  if  desired  we  can 
modify  £  uj^  so  that  it  vanishes  at  a  fixed  distance  outside  T. 


3.   APPLICATIONS 

In  this  section,  we  will  briefly  discuss  the  application 
of  the  extension  theorem  to  capacitance  matrix  and  iterative 
substructuring  methods.   We  note  that  the  main  idea  behind  the 


proof  of  our  first  result  is  due  to  Astrakhantsev  [  1  ] . 

There  are  occasions  when  a  linear  system  of  equations  can 
be  imbedded  in  a  larger  system  which  is  easier  to  solve.   For 
the  expanded  system  an  efficient  preconditioner  might  be  avail- 
able e.g.  if  the  geometry  of  the  corresponding  region  is  rec- 
tangular and  fast  Poisson  solvers  can  be  used;  see  Proskurowski 
and  Widlund  [24]  and  the  references  therein. 


Here  we  will 
for  selfadjoint  e 
The  region  Q  is 
venient,  simple  r 
sect  and  the  two 
symmetric  bilinea 
elliptic  problem, 
the  problem  is  st 


only  consider  the  solution  of  Neumann  problems 
lliptic  problems   by  capacitance  matrix  m.ethods, 
imbedded  in  A,  a  rectangle  or  some  other  con- 
egion.   The  boundaries  of  9.    and  A  do  not  inter- 
regions  are  triangulated  as  in  section  2.   A 
r  form  aj^(u,v)  is  associated  with  the  original 

To  avoid  nonessential  details,  we  assume  that 
rictly  elliptic  in  the  sense  that 


We  also  assume  th 
similarly  well  be 
considerable  free 
boundary  of  A. 


<_  a^  (u,u)  . 

at  this  bilinear  form  can  be  extended   to   a 
haved  form   ayy(u,v).   We  note  that  we   have 
dom.  in  choosing  a  boundary  condition  on   the 


Th 
satisfy 
a  confo 
of  sect 
to  solv 
meters 
Neumann 
def init 
matrix 
denoted 


,m. 


e  space  H  '(A),  or  an  appropriate  subspace  of  functions 
ing  these  boundary  conditions,  is  approximated  by  V  (A), 
rming  finite  element  spacg  satisfying  the  conditions 
ion  2.   The  restriction  V    {Q)    of  this  space  to  9.    is  used 
e  the  Neumann  problem  on  Q .   We  note  that  all  the  para- 
associated  with  V  (C)  are  determined  when  the  discrete 

problem  is  solved.   Denote  by   A  the  resulting  positive 
e,  symjT.etric  stiffness  matrix  and  by  nj-j  the  order.   The 
corresponding  to  the  problem  on  the  larger  region  A  is 

by  B  and  its  order  is  N,  . 


One  can  think  of  the  capacitance  matrix  method  as  a  pre- 
conditioned conjugate  gradient  method  in  which  the  trivially 
expanded  system 


A   0 
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is  solved  using  B  as  a  preconditioner.   As  always  the  conver- 
gence of  the  method  is  determined  by  the  spectrum  of 
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X  B  (f) 


or  equivalently  by  the  stationary  values  of  the  Rayleigh  quotient 
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Such  a  bound  is 


The  Nj^^-n,  linear 


g  the  extension  theorem^ 
ssociated  with  A  \  Q  is,  by 
tion  of  those  associated 


The  extension  theorem  therefore  provides  a  bound  of  the 
strain  energy  on  A  in  terms  of  that  on  Q    and  the  proof  is  com- 
pleted.  We  note  that  numerical    experiments     reported   in 
Proskurowski   and   Widlund  [   ]   show   that   'his  bound,   and 
thus  the  constant  in  the  extension  theorem.,  .   -eriorates  if 
triangles  become  very  thin. 
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structuring  metho 
discussion  of  thi 
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dimensions  higher 
belong  to  H^,  we 
To  simplify  the  n 
the  original  regi 
To,  the  intersect 
aries  of  oQ-[_,  Q2 
respectively.  We 
boundary  of  Q. 


to  a  discussion  of  an  optimal  iterative  sub- 
d;  cf.  Bj2$rstad  and  Widlund  [  7]  for  a  detailed 
s  and  similar  methods.   We  begin  by  discussing 
elements,  extending  our  previous  results  to 
than  two.   Since  Lagrangian  elements  do  not 
confine  ourselves  to  second  order  systems.      - 
otations,  we  also  concentrate  on  the  case  where' 
on  n   is  the  union  of  tv/o  regions  P.]_  and  ^2^^^ 

and  Qn.   The  bound- 
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The  variational  formulation  is  then 
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VqCP) 
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(3.1) 


where  V-_{_n)  is  the  subspace  of  functions  in  v  (Q)  which  vanish 
on  r,  *-"  ^2*   ^^^  again  assume  that  a^(u,v)  is  uniformly  elliptic 
and  selfaajoint. 


tend  an  element  of  Vt,(n^,r. 

:o  ah 


In  the  proof,  we  need  to    _ciiia  an  cxcmc=iiL.  wj.  v^  \jt.  ,1-1/, 
the  subspace  of  V^Cfi,)  of  functions  vanishing  on  r,,   to  an 
element  of  V}^(r2).   We  can  accomplish  this  by  extending  the  given 


0 


element  by  zero  in  CTl      and  then  use  the  construction  of  section  2 
to  obtain  values  in  P.^. 

The  linear  system  of  equations  corresponding  to  (3.1)  has 
the  form 
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where  K  is  positive  definite  and  symmetric-   The  matrix  K, , 
represents  the  couplings  between  pairs  of  degrees   of   freedom 
in  9.-^  ,     K-j^o  couplings  between  pairs  belonging  to  Q^  and  r3 
respectively   etc. 

We  make  the  assumptions  that  the  discrete  problems  on  the 

subregions,  with  some  appropriate  boundary'  condition  added  on 

r3,  can  be  solved  exactly.   We  can  therefore  set  b,  and  bn  equal 

to  zero.   By  block  Gaussian  elimination,  the  problem  (3.2)  is 
reduced  to 
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X3  =  b3 


(3.3: 


This  system,  is  preconditioned  by  S  ^  '  ,    a  matrix  related  to  a 
problem  on  fi-,  as  follows.   By  using  a  natural  boundary  condition 
on  r3  as  the  only  nonhomogeneous  data  and  a  zero  Dirichlet  con- 
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where  the  elements  of  K 
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(3.4: 


n  (*H  '*-i)  '  with  (}).  ,4  . 
"1   1   J  ID 


have  the  form  a 
basis  functions  associated  with  degrees  of  freedom  on  r  . 

The  problem  reduces  to 
(1 


5(1^x3  =  (K33 


T   -1 
-  K^3K^^K^3)X3  = 


Since  S  =  S 


(1) 


+  S 


(2) 


where  S 


(2) 


as  S^-^),  it  is  not  surprising  that   S^^) 
conditioner  of  S. 


is  constructed  in  the  sa-mc  'v;ay 
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To  provide  upper  and  lower  bounds  for  the  associated  gener- 
alized eigenvalue  problem,  we  consider  the  Rayleigh  quotient 
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harmonic  and  we  similarly 
of  the  solution  corres- 
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e  V^(r^^,r  )   which  also  satisfies 
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Therefore 

^n2^V^h'  V^h^  =  ^fi^^^^'^h^  ■"   ^fi2^^h'^h)  '  ^^h^^o'^'2'' 

i.e.  u,  has  the  least  strain  energy  of  any  element  in  ^^(^2/^2) 

with  the  given  boundary  values  on  T,.  In  view  of  this   the 

extension  theorem,  as  modified  in  this  section,  can  be  used  to 
complete  the  argument. 

We  conclude  with  some  comments  on  a  non-Lagrangian  case. 
To  be  specific,  we  consider  the  Dirichlet  problem  for  the  bi- 
harmonic  equation  approximated  by  the  18  degrees  of  freedom, 
reduced  quintic  element  discussed  in  Ciarlet  [13].   The  inter- 
polation conditions  are  associated  with  the  values  of  the  func- 
tion and  all  its  first  and  second  derivatives  at  the  vertices 
of  the  triangles.   Here  we  principally  v;ant  to  make  the  point 
that  there  are  two  different  ways  of  obtaining  an  H^(n)  exten- 
sion from  the  subregion  Pc-,       to  9^2      ^^"^  therefore  two  domain 
decomposition   algorithms.   For   simplicity,  we  assume  that  T^ 
is  straight. 

In  the  first  variant  all  the  degrees  of  freedom  associated 
with  r-,  are  shared  by  the  finite  element  functions  on  Q-^   and  ^2* 
The  proof  of  the  optimality  of  the  dom.ain  decomposition  algo-    ■; 
rithm  carries  over  directly  from  the  case  discussed  above.   In   ^ 
the  iteration,  approximate  values  of  six  parameters  per  vertex 
on  r3  are  obtained.   In  the  second  variant,  we  keep  the  values 

of    ^   free  when  extending  the  discrete  biharmonic  function 

9n^ 
from  Q,  to  P.^.   Since  this  corresponds  to  the  removal  of  con- 
straints we   have  no  further  problem  establishing  the  necessary 
bound.   From  an  algorithmic  point  of  view   there  is  a  possible 
benefit.   A  detailed  examination  of  the  linear  algebra  involved 
reveals  that  we  only  obtain  five  nontrivial  residuals  per  vertex 
on  r^,  corresponding  to  the  parameters  shared  between  fi^  and  ^2' 
We  note  that  the  values  of  d'^u,/dn'^      at  the   vertices   do   not 
affect  the  values  of  u^  and  its  gradient  on  F-j.   A  solution 
obtained  by  this  second  variant  will  thus  belong  to  H^  but  in 
general  it  will  differ  from  the  standard  discrete  solution 
obtained  by  applying  the  finite  element  method  on  the  entire 
region  Q. 
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